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The magnetism by the edge states in graphene is investigated theoretically. An instability 
of the pseudo-spin order of the edge states induces ferrimagnetic order in the presence of the 
Coulomb interaction. Although the next nearest-neighbor hopping can stabilize the pseudo-spin 
order, a strong Coulomb interaction makes the pseudo-spin unpolarized and real spin polarized. 
The magnetism of the edge states makes two peaks of the density of states in the conduction 
and valence energy bands near the Fermi point. Using a continuous model of the Weyl equation, 
we show that the edge-induced gauge field and the spin dependent mass terms are keys to make 
the magnetism of the edge states. A relationship between the magnetism of the edge states and 
the parity anomaly is discussed. 



CD 



o 



> 
o 



o 

00 

o 



X 



1. introduction 

The electronic properties of graphene^""^ have at- 
tracted much attention mainly because of its relativistic 
character of low energy electronic excitation. The energy 
band structure of graphene consists of two Dirac cones at 
the K and K' points in the fc-space. The electron dynam- 
ics around each Dirac point is approximated by the Weyl 
equation which describes a massless particle. When we 
consider solutions of the Weyl equation for a finite (or 
semi-finite) graphene cluster (ribbon) with the zigzag 
edges, the spatially-localized edge states'^ exist around 
the Fermi energy.^ The existence of the edge states de- 
pends on the shape of the edge for a graphene cluster. For 
example, the zigzag edge yields the edge states while the 
armchair edge does not. The energy dispersion for the 
edge states which appears only between the two Dirac 
points smoothly connects to the energy dispersion of 
the delocalized states.^ A large local density of states 
(LDOS) by the edge states may induce magnetism* and 
superconductivity^ near the zigzag edge. The existence 
of a mass for the Dirac particle is an important issue be- 
cause it gives rise to an energy gap at the Dirac points 
and relates to the ordered states. We examine a mecha- 
nism that the Coulomb interaction makes a mass and a 
magnetism. 

Fujita et al. discussed the effect of the Coulomb inter- 
action on the edge state, in which the electronic spins 
are localized at the edge to form a ferromagnetic state 
at one zigzag edge and another ferromagnetic state with 
the opposite spin at another zigzag edge.* The occur- 
rence of the magnetism is investigated by first princi- 
ples calculations, too.®"^" If the ferromagnetic state ap- 
peared at one edge, we would expect that two peaks for 
up and down spin states appear in scanning tunneling 
spectroscopy (STS). However this situation seems to be 
inconsistent with the STS measurements in which they 
observed only one LDOS peak near the zigzag edge of 
graphite. ^*"^^ Klusek et al}^ found a peak of LDOS in 
the energy range of 20-250 meV above the Fermi level at 
the edges of circular pits on graphite surface. Kobayashi 
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et al}^'^^ and Niimi et a/.^*'^'^ independently observed 
a peak in the LDOS below the Fermi energy by 20 ^ 30 
meV. Since the peak appears only conduction (Klusek et 
al.) or valence energy band (Kobayashi et al. and Niimi 
et al), it suggests that the edge states do not make a 
magnetism. Thus it is an interesting problem for under- 
standing the occurrence of the ferromagnetic order at the 
edge in the presence of the Coulomb interaction. When 
we see the calculated results of Fujita et al, the polar- 
ized spin appears for the edge states at a much small 
on-site Coulomb interaction U value compared with the 
nearest neighbor interaction 70 (see Fig. 5 in Ref. 4). 
Although we reproduce their results numerically, the re- 
sults are very surprising. A possible reason why the spin 
ordering occurs for such a small U is due to a special fact 
that the wavefunction of the edge states has an ampli- 
tude only one of the two sublattices (A and B) and thus 
the nearest neighbor interaction is suppressed. When we 
introduce the next nearest neighbor interaction, 7„, the 
spin polarization around U /70 ^ disappears and spin 
ordering appears from finite values of U depending on 
7„, which we will show in this paper. 

The wavefunction for two sublattice structure is re- 
ferred to as the "pseudo-spin" . An edge state can be de- 
scribed by a pseudo-spin polarized state. ^ A pseudo-spin 
structure gives a rich variety of interesting physical phe- 
nomena not only the edge states but also the extended 
states. For example, the absence of backward scattering 
mechanism is relevant to this pseudo-spin nature,*^' 
in which a 27r rotation of a pseudo-spin wavefunction 
around the K-point in the two-dimensional Brillouin zone 
does not gives the original wavefunction but gives minus 
sign to the wavefunction. Thus, the pseudo-spin is quite 
similar to the real spin in the real space. In this paper, we 
show that the pseudo-spin also plays an important role 
for the magnetism (or real spin) of the edge states, which 
is shown by a numerical analysis of the lattice model and 
by a analytical study of the Weyl equation. 

This paper is organized as follows. In Sec. 2 we explain 
the model Hamiltonian and introduce symmetric and an- 
tisymmetric variables for the pseudo-spin. In Sec. 3 we 
show numerical results for the ground state of the model. 
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In Sec. 4 we use a continuous model to examine the mech- 
anism of the magnetism of the edge states. Discussion 
and summary arc given in Sec. 5. 

2. pseudo-spin representation of Hamiltonian 

The Hamiltonian is given by Hoi + + Hu where 
Hos = -7oJ2{i,j)(^l,i(^s,j =1,1) is the nearest- 
neighbor tight-binding Hamiltonian (70 ~ 3cV is the 
hopping integral), and Hu is the Hubbard on-site in- 
teraction. Hu is written as 



Uu = U ^n|(r)n4(r), 



(1) 



where U is the on-site energy and n-|-(r) (n|(r)) is the 
density operator of up (down) spin electron at site r. 
Since the hexagonal lattice consists of two sublattice, A 
and B, Tiu is given as a summation over unit cells as 



Hu = UY^ n|,p(r„)n4,p(r„), 



(2) 



=A,B 



where n-f^p(r„) (n|^p(r„)) is the density operator of up 
(down) spin electron at p-sublattice (p = A, B), and r„ 
denotes the position of a unit cell. For a unit cell, we in- 
troduce a density and a magnetization at r„ as np(r„) = 
nT,p{ru) + ni^p{ru) and mp(r„) = nt,p(r„) - n4,p(r„), re- 
spectively. Hereafter r„ for each variable is omitted for 
simplicity. From rip and mp, we define density and mag- 
netization for a unit cell as 



n = nA + nB, m = mA + "^B• 



(3) 



n and m are symmetric with respect to the sublattice. 
Here we introduce pseudospin order and antiferromag- 
netic order for a unit cell, 



Pn 



riA - riB, 



niA - thb, 



(4) 



which are anti-symmetric with respect to the sublattice. 
Pn (Pm) represents charge (spin) polarization within the 
hexagonal unit cell. T-Cu can be rewritten in terms of 

n, m,Pn and Pm as 

This representation of Hu shows that not only non- 
vanishing magnetization ((m) 7^ 0) but also antiferro- 
magnetic ((m) = and (pm) ^ 0) or ferrimagnetic 
((m) ^ and (pm) 7^ 0) spin configuration are favored 
to decrease Hu where (O) denotes the expectation value 
of operator O for the ground state. 

By applying the mean-field approximation to Eq. (2), 
'Hu = UJ2^^ p{ni,p)ni^p + nt,p(nj.^p) - {ni,p){ni,p), the 
Hamiltonians for up and down spin electrons in graphene 
are given as follows: 



Hf 



Ho\ + 



-E 
2 ^ 



{UA - JTIa) 



(nA + niA) 




\ (n^.A 
(riB - tob)/ Vn|,B 





(riB + tob) 



(6) 



The Hamiltonians of Eq. (6) can also be rewritten using 
Eqs. (3) and (4) as 



Hi=Hoi + -J2 [{{Pn) - iPm)) a. + ((n) - (m)) /] 



"-T,A 

'^i.A 



(7) 

where <Tz = diag(l, —1) and I = diag(l, 1). The pseudo- 
spin variables arc proportional to ctz and affect magneti- 
zation of the edge states as we will show in the following 
sections. 



Table I. Parities with respect to spin and pseudo-spin 

Pn Pm n m (Tz CoupUng 



Spin (T-^i) + - 

Pseudo-Spin (A<-+B) — - 



In Table. I, we show the parities of j^^j,, Pm,n,m and 
for changing the direction of spin and pseudo-spin. The 

interaction terms in Eq. (7) are invariant with respect 



to pseudo-spin parity: p„ 



Pn'; Pn 



m ^ m and az 
spin parity: p„ — 



—(Tz, and we have W| Hi for 



Pn , Pm 



-Pn 



-m and 



— *■ • Since Hu is even parity with respect to spin and 
pseudo-spin, it is expected that a ground state is realized 
by spontaneous symmetry breaking if Hq is symmetric, 
too. Ho is symmetric with respect to spin but asymmet- 
ric with respect to pseudo-spin due to the presence of 
the zigzag boundary. This can be explained as follows. A 
magnetic field, B*"", breaks the spin degeneracy of the 
ground state by the Zeeman term, — B**™ • m, which is 
odd parity with respect to spin. Similarly, we can de- 
fine a pseudo- magnetic field, H^, that couples to az (see 
Eq. (10)) and breaks the degeneracy of the pseudo-spin 
parity. In fact, it can be shown that B'l appears at the 
zigzag edge^ so that Ho can induce a pseudo-spin or- 
der, (pn) ^ 0, near the edge. In the next section, we will 
show numerically that Ho of a zigzag nanotube breaks 
the pseudo-spin parity of Pn, which is an important to 
obtain magnetism for the total Hamiltonian, H^ +Hi. 

3. numerical results and analysis 

In this section, we show numerical results for the 
ground state of H^ + Hi. {n/4), (pn), {Pm), and (m) are 
plotted for (50, 0) zigzag nanotube with length L « 4nm, 
and LDOS curves are calculated for (100,0) zigzag nan- 
otube with length L « 20nm. We set the origin of the 
Fermi energy Ey = as Hs = U /2. 

In Fig. 1(a), we plot (n/4) and (p„) in the case of 
[/ = 0. Since ?7 = 0, (m) = (p,„) = 0. The solid 
(dashed) curves are the results for the Fermi energy 
= -0.01 eV (-hO.Ol eV). (n/4) and (p„) are mod- 
ulated near the edges and their difference from the con- 
stant values ((n/4) = 0.5 and (p„) = 0) is due to the 
presence of the edge states.^ The wavefunction of the 
edge states is localized near the edges so that (n/4) 
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is different from 0.5 (i.e., half filling n = 2) only near 
the edges. Moreover, the wavefunction of the edge states 
is polarized about the pseudo-spin and (p„) is nonzero, 
too. (n/A) and are unstable against a small change 
of Ep due to the flat energy band of the edge states. 
In Fig. 1(b), we plot the LDOS at L w 0.5 nm in the 
case of = 0. The spin up edge states and spin down 
edge states are degenerate in the case of [/ = so that 
they make a sharp LDOS peak at Ep = 0. To show the 
(smooth) LDOS curve, we put a constant width (0.05 
eV) for each state. 



(a) 



Ef = ±0.01 eV 



(a) 



(b) 



U = 



Ey = +0.01 eV 
Ey = -0.01 eV 




1 1 

E (eV) 



Fig. 1. (a) Self-consistent solution of {n/4} and (p„} for (7 = 
and Ey = ±0.01 eV. Due to the presence of the edge states 
consisting a fiat energy band at E-p = 0, {n/4) and {p„) are 
unstable against the small change of Ep. {0{ru)) depends only 
on the distance from an edge due to rotational symmetry around 
the axis of the tube, (b) Corresponding LDOS curve at L si 0.5 
nm. The degenerate spin up and down edge states make a sharp 
peak at Ep = in the case oi U = 0. 



In Fig. 2(a), we plot self-consistent solution of (n/4), 
(m), (pn), and (pm) in the case of [/ = 79. The result 
shows that, (n/4) « 0.5 and (p„) « hold at each hexag- 
onal unit cell, and (m) and (pm) become nonzero near the 
zigzag edges. Since Hu can stabilize (n) and (p„) accord- 
ing to Eq. (5), the pseudo-spin polarization {{pn) ^ 0) 
which exists for U = disappears for the ground state 
due to a finite value of U. The corresponding LDOS curve 
is shown in Fig. 2(b). Because of U, the spin up (down) 
edge states are shifted above (below) the Fermi energy so 
that the ferrimagnetic order ((m) ^ and (pm) 7^ 0) ap- 
pears and that two LDOS peaks appear around Ep ~ 0. 
The LDOS curve in the case oi U = is also shown in 
Fig. 2(b) for comparison. 

The numerical results can be explained qualitatively 
using Eq. (7) as follows. When (p„) = and (pm) > 0, 
the energy of a spin up electron is shifted below the Fermi 
energy for a pseudo-spin up state (ctz = 1)- On the other 
hand, the energy of a spin down electron is shifted above 
the Fermi energy for the same state (a pseudo-spin up 
state). Thus, when Ep = 0, the ground state has a finite 
positive value of (pm) and (m), which lowers the energy 
of a spin up electron due to the last term in Eq. (7) and 
stabilizes the ground state configuration further. In fact, 
the second terms of the right-hand side of Eq. (7) give 
rise to an energy gap in the energy spectrum. As we will 
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Fig. 2. (color online) (a) Self-consistent solution of (n/4), (m), 
(pn) and (pm) for U = -fo and Ep = iO.OleV. The spin around 
the edge is polarized and |(pm)[ ~ l(™)l shows that spin is a 
ferrimagnetic configuration, (b) The LDOS curve at L fx 0.5 
nm shows that spin up edge states (denoted by up-arrow) and 
spin down edge states (denoted by down-arrow) make two sharp 
peaks (solid curve) near Ep = when U ^ 0. 



show using a continuous model for Ti.Qs in Sec. 4, the 
appearance of the gap will become more clear since the 
term proportional to CTz acts as a mass term of Dirac 
fermion. It is noted that (p„) = is consistent with the 
presence of a gap, and non- vanishing (pm) gives different 
signs of the mass terms for spin up and down electrons. 

Next, we consider the next-nearest neighbor (nnn) 
hopping, 7„, which is an intrinsic perturbation to the 
edge states. In the previous paper, we showed that the 
nnn interaction, Tinnn, gives a finite energy bandwidth 
to the edge states, W = 7„.^^ 7„ « O.I70 is obtained 
by first-principles calculation using the local density ap- 
proximation.^^ A finite energy band width of the edge 
states suppresses the above mentioned (pn)'s instabihty 
with respect to a small change of Ep (see Fig. 1). 

To see the relationship between (p„) and (pm) in detail, 
wc first define the net pseudo-spin order P„ and the aver- 
aged antiferromagnetic order Pm- Pn = J2u l(Pn(''"))l/'^ 
and Pm = X]«(Pni(''"))/'^ where the summation is taken 
over all hexagonal unit cells. Since not only the edge 
states but also extended states can contribute to Pi 
{i = n,m), we consider the difference between Pi for 
a tube with the zigzag edges (P/"*^") and that for a cor- 
responding periodic torus system (p^'o^us^ which does not 
have edge. P^^s° = P^^'' - P^°™" can be used to show 
the magnetism for the edge states. 

In Fig. 3(a), we plot P^^"" (solid curve) and Pf^'' 
(dashed curve) as a function of C//70 for 7„ = (black), 
O.I70 (red), and O.270 (blue). When 7„ = 0, no pseudo- 
spin order P^'^sc ~ q fgj. g^j^y positive value of U ^ while the 
antiferromagnetic order P^yf''^'^ increases until f//7o ~ 2.1. 
P^'^se decreases when C//70 > 2.1. However Fujita et al. 
showed that the magnetism due to the extended states 
increases.^ In case of a finite value of 7„, the antiferro- 
magnetic order is suppressed P^^° ~ up to a finite 
value of U and the magnetism (in the case of 7„ = 0) 
appears discontinuously above the critical value of U . On 
the other hand, the pseudo-spin order appears P^'^s'' ^ 
below the critical value of U . Thus, when Ep = 0, we 
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Fig. 3. (color-online) (a) Self-consistent solution of Pm^" (solid 
curve) and P^''^" (dashed curve) as a function oiU/-fo for 7n = 
(black), O.I70 (red), and O.270 (blue) in the case of Ep = 0. 7,j 
controls the appearance of the pseudo-spin order {Pn'^^'^) and 
of the ferrimagnetic order (P„ '^°). (b) The LDOS curves for 
7n = O.I70 (solid curve) and 7„ = O.270 (dashed curve) when 
(7 = 70 and Bp = 0. (c) Pm^" (solid curve) and P^'^^" (dashed 
curve) for 7^ = O.270 with Ep = — O.I70. This corresponds to 
that Ep is located at the center of the edge energy band. The 
small step on Pm^^ (at C//70 = 1) is due to a finite diameter of 
a tube, (d) The corresponding LDOS curve when f7 = 70 in (c). 
The peak structure is sensitive to the position of Ep . 



see that 7„ controls the occurrence of the pseudo-spin 
order (-P°'^^°) and the antiferromagnetic order {P°^^'^) 
exclusively. In Fig. 3(b), we plot the LDOS curves when 
i7 = 7o for 7„ = O.I70 (solid curve) and 7„ = O.270 
(dashed curve). When 7„ — O.I70, magnetic order is re- 
alized so that there are two peaks in the LDOS curve. 
When 7„ = O.270, the pseudo-spin order P°'^s° is real- 
ized and there is one peak below the Fermi energy in the 
LDOS curve. In the case of J7 = 0, the peak position 
appears a,t E — — 7„.^* The pseudo-spin order shifts the 
peak position above E = — 7„ due to (J7/4)^j. {pn)crz 
in Eq. (7). This is a possible reason why Kobayashi et 
al}^' and Niimi et al}'^' observed a peak in the LDOS 
below the Fermi energy by 20 ~ 30 meV not by 7„ « 0.3 
eV.^^ Since the most localized edge states have p„ = ±1 
at a unit cell of the edge site and the energy is given 
by -7n, we see that Hnnn = -Jn\pn\- Thus, the p„ de- 
pendent energy density at a unit cell of the edge site, 
Uu + Unnn = (J7/8)p^ - 7„ |p„ | -|- . . ., may become a neg- 



ative value when the ground state shows a pseudo-spin 
order, (p„) ^ 0. In fact, in the case of 7„ = O.270, the 
magnetization disappears even for [/ « 79. 

In Fig. 3(c), we plot P^^° (soHd curve) and Pf^"" 
(dashed curve) as a function of C//70 for 7„ — O.270 with 
E-p = — O.I70. The results show that, the critical U value 
decreases as compared with the case of E-p ~ shown by 
the blue curves in Fig. 3(a). The corresponding LDOS 
curve for [/ = 70 is plotted in Fig. 3(d). We see that 
there are one sharp peak and a broaden peak. Although 
P^'^sc p^tig® when C/ = 70, this two peaks structure is 
not so clear as the two peaks in the case of Ep = with 
7„ = O.I70. 

4. Continuous model 

To understand the edge magnetism in detail, we solve 
Eq. (7) analytically by means of a continuous model. It 
will be shown that the magnetic order is explained by a 
gauge field for the edge states and spin-dependent mass 
term. 

The low energy states around Ep — consist of elec- 
trons near the K-point and K '-point. Since the K-point 
and K '-point are related to each other by time- reversal 
symmetry, it is sufficient to consider only the K-point 
when B**™ = 0. Then, the low energy Hamiltonian is 
given by replacing Hqs in Eq- (7) with TYk as 

U 



where is given by 
Hk = 



VpCT ■ (p - 



(Pm))o-2 



(8) 



(9) 



In Eq. (8), Hs (s =Tji) operates on a two-component 
wavefunction, ipf = '(V'a.s' V'b s)' where i/'as ^'^^1 i/'bs 
are the pseudo-spin up and down states, respectively. In 
Eq. (9), vp is the Fermi velocity, p = (px^Py) is the mo- 
mentum operator, and <t = {ax,cry) where ai (i = x, y, z) 
are the Pauh matrices. A'i(r) = (A^(r), AJ](r)) is a field 
that is induced by a defect (edge structure) in the hexag- 
onal unit cell and is referred to as the deformation- 
induced gauge field. ^"'^^ The corresponding deformation- 
induced magnetic field, i?5(r) = dxA'^{r) — dyA^{v) , cou- 
ples to the pseudo-spin, cr^ . This is shown by TYk squared, 

Hk = 4 {(P + A'i(r))2 -f hB'^{v)az] , (10) 

where az selects the direction opposite to B'^{r) in or- 
der to decrease the energy. In the previous paper, we 
derived A'^{y) and B'^(y) for zigzag edges. ^ A'^(r) ap- 
pears when we cut the hexagonal lattice at dashed lines 
in Fig. 4. Kusakabe et al. discussed the edge state mag- 
netism for the two possible edge structures,^ that is, the 
zigzag edge and the Klein edge.^^'^^ In Fig. 4, we show 
the corresponding A''(r) and P^ (r) for a zigzag edge (a), 
a Klein edge (b), and a graphene cluster with a zigzag 
edge at one edge and a Klein edge at another edge. The 
direction of A'^(r) for the Klein edge is opposite to that 
of the zigzag edge. In all cases, we can explain the edge 
structure within the same frame. 
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o +$0/4 • -$0/4 



Fig. 4. Configuration of A'^(r) and -Bz(r) for (a) the zigzag edge 
and (b) the Klein edge. When we cut the graphene sheet at 
the dashed line, the edge and Ai(r) = (Ax{y),0) appears. The 
direction of A'5(r) for the Klein edge is opposite to that of the 
zigzag edge, (c) A graphene with the zigzag edge and the Klein 



In the following, we obtain the wavefunction and the 
energy eigenvalue E for the following Hamiltonian: 

(Wk +m«a.) V^^,,(r) = E^l^,{r), (11) 

where the mass term is given by 

rris = j((Pn> T (Pm))- (12) 

The mass term depends on real spin, that is, the neg- 
ative (positive) sign in front of (pm) is for s =t (i). 
Equation (11) describes Dirac fcrmion having a mass, 
nis: whore the dimension of nig is energy here. In ob- 
taining Eq. (11), we neglect if{U/4){rn) of Eq. (8). The 
neglected term docs not couple to the psc;udo-spin and 
only shifts the energy position of each state so that it is 
not important for our discussion. Further, we neglect the 
r dependence of in order to simplify the argument. 

Since there is translational symmetry along the edge, 
the eigenfunction of Eq. (11) can be expressed by 

Vp.,«(r) = A^exp (i^x) e'^^^) (^t^j^J) , (13) 

where s =T)i, a-ud x (y) is parallel (perpendicular) 
to the zigzag edge.^ The unknown functions, G(y) and 
gs{y), and E can be determined by putting Eq. (13) into 
Eq. (11). We obtain 

Px + My) + ^^9s{y) =Dcos'h{2g,{y) + /,), (14) 



n—G{y)=Dsmh{2g,{y) + f,), 



(15) 



where variables D and fs are defined respectively as 

D = ± — ^/E^ -ml (16) 



Vf 



tanh(/s) = - 



E 



(17) 



Here, we consider a localized wavefunction and put 
G{y) = \y\/^ into Eq. (15) where ^ is localization length 



of the edge state and the zigzag edge is located aty = 0. 
Then we get 



^9s{y) + fs 



sinh- ( — 



+ sinh 



AD, 



(y<o) 
{y > 0). 



(18) 



Next, we integrate Eq. (14) with respect to y from —^g 
to ^g. By considering ^ 0, only singular functions of 
A^{y) and Eq. (18) around y = survive, and we get 

?9 



— sinh 



AUy)dy- 



h 

Using Eqs. (18) and (19), we see from Eq. (14) that 

P: 



(19) 



D 



cosh A'^{y)dy 



(20) 



holds except very close to the edge. From Eqs. (19) and 
(20), we see that ^ in G{y) is given by 



I = -Px tanh 




AUy)dy 



(21) 



This result is surprising in the sense that ^ in the presence 
of U is identical to ^ for U = 0.^ The mass term would 
affect ^, but it is not the case. The reason for this will be 
discussed elsewhere. Finally, we get the energy eigenvalue 
from Eqs. (16) and (20), 

{ypPx)"^ 



E' = mi + 



cosh^ (^fj^^ AUy)dy) ' 



(22) 



The energy dispersion relation for the edge states of 
Eq. (22) is similar to the relativistic energy dispersion 
relation for the extended state, E^ = ml + {v-ppxf- 

In Eq. (17), we see that the sign of fs depends both 
on the signs of mg and E. To obtain a ground state, 
we first consider the valence states E < Q. Then, we 
have fs = sign(ms)|/s|. Using Eq. (19), we can rewrite 
Eq. (18) as 



9s{y) 



A'^{y')dy' - -sign{ms)\fs\ (y < 0) 



A';^{y')dy' - -sign{ms)\fs 



By putting Eq. (22) into Eq. (17), 

f_i AUy)dy 



ha- 



iy > 0). 

(23) 

e |/,s| « 



when 



0. Since 



III Al{y)dy 

J-\ A^{y)dy » for the zigzag edge (see Fig. 4(a)), 
we get from Eq. (23) 

Al{y')dy' (t/<0) 
{y > 0), 



9\{y) 



(24) 



when m| < 0. The wavefunction of this spin up state ap- 
pears near the edge consisting of A-atoms {y < 0) in the 
valence band {E < 0): ipp^^iiy < 0) oc exp(-|y|/^)'(l, 0). 
The wavefunction at the edge consisting of B-atoms 
{y > 0) is pseudo-spin unpolarized and the amplitude 
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is negligible due to the normalization of the wavefunc- 
tion. Similarly, for mi > 0, we have 



(b) (p„> = E 



9i{y) 







(y<o) 
{y > 0). 



(25) 



The corresponding wavefunction of spin down state is 
pseudo-spin down state appearing only near the edge 
consisted of B-atoms (y > 0). It is noted that the spin for 
a conduction edge state {E > 0) is opposite to that of a 
valence edge state. Thus, we obtain local ferrimagnetism 
near the zigzag edge. 

In Fig. 5, we show how the magnetism appears around 
the zigzag edge. In Fig. 5(a), if t/ = 0, the opposite di- 
rection of pseudo-spin (the edge state) appears both for 
y > (fT^ = —1) and y < (cFz = !)• However when 
U ^ 0, due to the mass term the edge state exists only 
one of the two sides in Fig. 5(a). When pseudo-spin or- 
der is suppressed {{pn) = 0) and antiferromagnetic order 
appears ((pm) 7^ 0), we get a situation of a different sign 
for m-f and mi {mi = —mi) from Eq. (12). In this case, 
up (down) spin edge state appears for y < (y > 0) for 
the valence band according to Eqs. (24) and (25), while 
down (up) spin edge state appears for y < (y > 0) for 
the conduction band. It is consistent with the numerical 
result of Fujita et al. for Ep = O,** in which the electrons 
are occupied only for the valence band. If Ep is shifted 
from Ep = then (pm) will disappear. 

When pseudo-spin order appears {{pn) 0) and an- 
tiferromagnetic order is suppressed {{pm) = 0), we get 
another situation that the sign of mi and mi are the 
same (m-f = m^) from Eq. (12). In this case, the energy 
levels for up and down spin edge states are degenerate 
even for J7 ^0 (see Fig. 5(c)). Thus the ground state 
is not spin polarized but pseudo-spin polarized. Even in 
case that mi ^ m^, spin up and down edge states both 
appear below the Fermi level so that the ground state is 
still spin impolarized if sign(TO|) = sign(m|). It is inter- 
esting to note that, in Fig. 5(c), the energy level position 
for y < can appear above (below) = when rrig > 
(?Tis < 0). The sign of ms, that is, ms < (m,, > 0) for 
y < and ms > {ms < 0) for y > 0, is consistent with 
the numerical results given in Sec. 3 

As we have shown in Sec. 3, the magnetism of the edge 
states is affected by the nnn hopping that can stabilize 
{Pn)- In the continuous model, we showed in the previous 
paper that the nnn perturbation works as 



= 7J^B^{r)a, 



(26) 



for the edge states where £ = 3acc/2.^ Wnnn is pro- 
portional to az so that Wnnn appears as an additional 
term for the mass. If 7„ is sufficiently large then we have 
TOj = mi and magnetism disappears. This is consistent 
with the numerical results given in Sec. 3. 

5. Discussion 

A magnetism of the edge states would give rise to two 
LDOS peaks since only spin up (or down) edge states are 
located below the Fermi energy to give a spin polariza- 
tion in the ground state. Although we have examined this 



mi 



az = l 

(Tz = l Oz = -1 (C) (pm) = £ 



0^-mi 

<Tz = -1 



'Zigzag y 



A' 


B 


TOs < 


mg > 






az = l 


az = -1 



Fig. 5. (a) Deformation-induced gauge (magnetic) field produces 
asymmetry of pseudo-spin. The pseudo-spin polarized localized 
state (i.e., the edge states) appears in pair with respect to j/ = 0. 
(b) Due to the mass term, this symmetric character with respect 
to ?/ = is lost and the edge state can appear independently 
at J/ > and y < (spatial parity with respect to y — > —y is 
broken by the mass term) . When the sign of the mass is opposite 
with respect to spin, namely when m| = —mi, spin becomes 
asymmetric with respect to E = 0, and a local ferrimagnetism 
appears, (c) In case that = , spin up and down states are 
degenerate and both appear below the Fermi level when ms < 
(ms > 0) for 2/ < (2/ > 0) so that the ground state is spin 
unpolarized. 



mechanism using the Hubbard model, the appearance of 
two peaks seems to be a model independent consequence 
of the magnetism of the edge states. The LDOS near the 
zigzag edge of graphite has been measured by STS,^^~^^ 
but no experimental group has observed the two peaks 
in the STS data. It is possible that the position of the 
Fermi energy in these experiments is not suitable for the 
occurrence of the magnetism (see Fig. 3). Thus if we 
change of the Fermi energy, LDOS will give a split of the 
peak, which is an evidence that the edge states form a 
magnetism. 

We explained the magnetism of the edge states in 
terms of the spin dependent mass terms and the defor- 
mation induced gauge field. It is known that the mass 
and a gauge field in the Weyl equation induce the par- 
ity anomaly in the ground state. The mass term in 
Eq. (11) changes its sign under spatial parity with re- 
spect to y — » — y and V'^,s ~* ^xi'^^^s- The mass term 
violates the spatial parity and can induce a quantum 
anomaly in the ground state, which is referred to as the 
parity anomaly. By applying the formula of the parity 
anomaly^^ to our case, we obtain 



, , , , 1 B^{r) . , , 



(27) 



where $0 = 27rfi, is the flux quantum and correction may 
arise due to higher order derivatives of -B5(r). In the 
case of mi = —mi, we have magnetism, i.e., (pt(r) — 
Pi(i-)) + 0. Moreover, using /^,^;^ ^^^^B}Xr)(e-v = ±$o/4 
that will be derived in the following, the magnetization 
at the edge is estimated by (m) = ±1/4, which is good 
agreement with our numerical result, (m)o.L ~ ±0.25, 
shown in Fig. 2. Thus, we think that there is a close 
relationship between the edge states magnetism and the 
parity anomaly. In fact, the anomaly survives even in the 
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massless limit nis 0, which is consistent to the fact 
that an infinitesimal value of U gives a finite magnetism 
if we do not consider the nnn interaction (sec P^'^se 
Fig. 3(a) and Fig. 5 in Ref. 4). A graphene with the zigzag 
edge and the Klein edge can be used to know that the 
deformation induced magnetic field in a unit cell at one 
side of the edge is given by J^^.^ ^^^^B'^{r)d'^r = ±$o/4 
(see Fig. 4(c)). In this case, = J B'^{r)cfr is nonzero 
and the index theorem^'^ can be used to know B^{r). The 
thcorcni states that Hk possesses |$'^/$ol zero energy 
edge states. Since it is known that the number of the zero 
energy states in a (n, 0) nanotube is given by n,^^ then we 
can know that the flux in a unit cell at the zigzag edge, 
is given by ±$o/4. Here, we used 2x (2n$„) = ±n$o- 
The factor 2 comes from the time-reversal symmetry (the 
K and K' points) and the factor 2n is the total number of 
edge sites at the zigzag edge and the Klein edge. = 
$o/4 is consistent with the numerical result by Nakada 
et al?^ who demonstrate that an edge shape with three 
or four zigzag sites per sequence is sufficient to show an 
edge state. 

In summary, we have shown that the instability of the 
pseudo-spin order of the edge states induces ferrimag- 
netic order in the presence of the Coulomb interaction. 
The nnn hopping can stabilize the pseudo-spin order, but 
a larger value of U makes the pseudo-spin order unpo- 
larized and gives rise to a ferrimagnetic order. The ferri- 
magnetic order is sensitive to the Fermi energy position. 
In case that the pseudo-spin order is realized one peak 
appears in the LDOS near the zigzag edge, which is con- 
sistent to the existing experimental results. Using a con- 
tinuous model of the Wcyl equation, we showed that the 
deformation-induced gauge field gives rise to the mag- 
netism of the edge states if the mass terms have different 
sign for different spin edge states. 
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